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^ Abstract 
Oh 

We construct a supersymmetric generalization of the holographic dual of a fractional 
topological insulator found in [Ij. This is accomplished by introducing a nontrivial gauge 
field on the world volume of the probe D7 brane. The BPS equations are derived from 
,— I the K-symmetry transformation of the probe brane. The BPS equations are shown to 

reduce to two first oder nonlinear partial differential equations. Solutions of the BPS 
equations correspond to a probe brane configuration which preserves four of the thirty- 
Q two supersymmetries of the AdS^ x background. Solutions of the BPS equations which 

correspond to a holographic fractional topological insulator are obtained numerically. 
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1 Introduction 



The AdS/CFT correspondence, and more generally gauge/gravity dualities, provide a powerful 
tool for studying strongly coupled field theories in states with finite density, and hence might 
be useful in condensed matter physics (see [2], [3l H] and references therein). For example, 
superfluids [S], El E] and Non-Fermi liquids |S] can be realized in bottom-up models within 
this framework. Note that such bottom-up models are generically formulated in terms of four 
or five dimensional gravity theories coupled to scalars and gauge fields. It might or might 
not be possible to embed such models into string theory. Consequently, the dual field theory 
formulation in terms of elementary fields and a Lagrangian is in general not known or might 
even not exist. String theory embeddings of these bottom-up models can be realized by using 
probe branes in the background of D3-branes. In this case the dual field theory is known 
explicitly. For example p-wave superfluids [Ulinilll], Fermi surfaces jT2] and the gravity dual 
of a Quantum Hall Plateau transition [13J were investigated. 

In this paper we consider the dual gravity description of a fractional topological insulator. 
A (conventional) topological insulator is a proposed new type of quantum matter which is 
not adiabatically connected to an ordinary insulator. They are characterized by fully gapped 
excitations in the bulk and gapless boundary modes, whose vanishing mass is protected by a 
discrete symmetry. Topological insulators have been a very active field of research in the past 
couple of years (see [HI [151 llHl [Ej for reviews with references to the original literature). 

A very simple model for a Z2 time reversal invariant topological insulator ^2] is given by 
a Dirac fermion in 3-1-1 dimensions with a spatially varying mass 



One considers an interface across which m{x) jumps from a positive to a negative real value. 
For such an interface there is are massless localized fermionic degrees of freedom [121 1^ . 
realizing the gapless boundary mode of the topological insulator. The interface separates 
topological trivial and nontrivial phases of a Z2 topological invariant of the electron system 
(such as the Q term ^ E • B where = corresponds to the trivial Z2 insulator while 9 = ti 
corresponds to the Z2 non-trivial insulator). 

By analogy to the relation of the integer and fractional quantum hall effect we can ask 
whether it is possible to construct a time reversal invariant fractional topological insulator for 
which is a non-integer fractional multiple of vr. In three spatial dimensions such a realization 
of a fractional topological insulator is known [21] (see also [221 123]). The idea is that the charge 
carriers, i.e. electrons, are made up of partons which carry fractional charge. To ensure that 
the partons confine into electrons outside the topological insulator we have to add a statistical 
gauge field which is deconfined inside a topological insulator. 
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The structure of this paper is as follows. In section |2] we review the probe brane setup for 
a fractional topological insulator following [T] . Such a system breaks all supersymmetries due 
to the position dependent mass. In section [3] we present a field theory analysis that identifies a 
counterterm which restores some of the broken supersymmetries of the topological insulator. It 
is shown that the introduction of such a counterterm can be engineered by turning on a specific 
gauge field on the D7 brane probe, emulating ideas which were used in the construction of 
supersymmetric interface theories [Ml ESI EZl • In section |4] we construct the BPS equations 
from the /t-symmetry transformations of the probe D7 brane. The conditions that four of 
the thirty-two supersymmetries of the AdS^ x vacuum are preserved reduce to two first 
order nonlinear partial differential equation^ In section [s] we obtain solutions using various 
techniques. The solution which has the right asymptotic behavior to describe a topological 
insulator and is everywhere regular is obtained numerically. We discuss possible applications 
and generalizations of our results in section |6j Conventions, review material and technical 
details of the calculations are relegated to various appendices. 



2 A holographic fractional topological insulator 

In ^21] (see also [221 123]) the idea of a fractional topological insulator was introduced. In 
this model the charge carriers fractionalize into partons with fractional charge. Consequently, 
while the system is still time reversal invariant, the topological invariant does not have to be 
Z2 valued any more. Additional degrees of freedom (i.e. a strongly interacting gauge field 
coupled to the partons) ensure that outside the topological insulator the partons are confined 
and only integer charged electrons appear. 

A holographic realization of a fractional topological insulator was presented in [T], fIE\ . 
where the partons are realized by adding matter to a A/" = 4 SYM "phonon" bath. We give 
a brief review of this construction following [29] . The starting point is the intersection of D3 
and D7 branes in ten dimensional flat spacetime (with coordinates Xq,Xi, . . . ,Xg). 
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Table 1: D3/D7 system. The bullets denote Neumann boundary conditions. The D7 brane 
can be separated from the D3 in the 8-9 directions. 



^Note that a supersymmetric probe D7 brane preserves at most sixteen of the thirty two supersymmetries 
of the vacuum. 
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Whe choose the Nc coincident D3 branes to be located at X4 = ■ ■ ■ = X-j — Xg = Xg = 0. 
The massless modes in the D3-D3 string sector give rise to an A/" = 4 supersymmetric Yang- 
Mills theory with gauge group SU{Nc). From the point of the fractional topological insulator, 
these supersymmetric gauge theory degrees of freedom should be viewed as "phonons". 

Moreover, we also add fundamental matter degrees of freedom to our setup by embedding 
Nf coincident D7 branes. The D7 branes are aligned along the Xo,Xi,...,X7 direction. 
The world volume coordinates Xa {a e {0, 1, ... , 7}) can therefore be identified with Xa, 
i.e. Xa — Xa- If the D7 branes are located at Xg = Xg = 0, the lowest energy excitations 
of the D3-D7 strings describe massless matter. The symmetries of this configuration are 
5*0(4) = SU{2)r X SU{2)l rotations in the X4-X7 space and S0{2) = U{1)r rotations in 
the Xg — Xg plane. The massless modes in the D3-D7 string sector can be written in terms 
M = 2 hypermultiplets where the SU {2)r x U {1)r rotations correspond to the superconformal 
R-symmetry. Moving the location of D7 brane away from Xg = Xg = breaks the U{\)r 
symmetry and introduces mass terms for the hypermultiplets. 

In the large Nc limit, for A = NcQym ^ ^^e Nc coincident D3 branes can be replaced 
by the AdS^ x geometry, which is conveniently expressed as follows 

ds^ = ^ ^2^ i-dt^ + dxl) + (dxl + dxl + dxl + dx] + dXl + dX|), (2.2) 

where 



r 



2 - + xl + + a;2, p2 ^ ^2 ^ ^2^ ^2.3) 



In these coordinates the AdS^ boundary is located at r — > cxo. i? is the curvature radius of 
AdS^ and is given by i?^ = A From now on we will set i? = 1 and therefore = A. 

Let us now embed the Nf coincident D7 branes. In the present paper we will only consider 
the the probe limit (i.e Nf <^ Nc), where the back-reaction of the D7 brane can be neglected. 
Then the flavor degrees of freedom are described by embedding the D7 brane world volume 
along an (asymptotically) AdS^ x subspace inside the AdS^ x S^, with world volume 
coordinates Xa- 

Without loss of generahty, let us single out the field theory direction X3 and consider a non- 
trivial profile m{xs) for the mass of the hypermultiplets. On the gravity side, this corresponds 
to a non-trivial embedding of the D7 branes in the transverse (Xg, Xg) space of the form 

X8 = Xg(r,a;3), Xg = 0. (2.4) 

Since Xg 7^ 0, the U{1)r symmetry is broken. In principle, Xg can be any function of Xa with 
a e {0, . . . , 7}. In order to preserve the SO{A) = SU{2)l x SU{2)r symmetry Xg depends 
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only on r and not on Xa with a G {4, 5, 6, 7} explicitly. The mass m{x3) can be read off from 
the boundary behavior of X8(r, X3) 

lim Xs{r, X3) = m{xs) + O ( \] . (2.5) 

Besides the scalar functions Xg 9 there are also gauge fields living on the D7 branes. In case of 
Nf coincident D-branes the gauge fields are valued in U{Nf). The dynamics of Nf coincident 
D7 branes is given in terms of a Dirac-Born-Infeld (DEI) term, as well as a Wess-Zumino 
(WZ) term. As we will be interested only in the U(l) part of the U{Nf) world-volume gauge 
fields and scalars, the relevant part of their action reads 

Sd7 = Sdbi + Swz (2.6) 

with 

Sdbi = -Nffx-r j d^'x^- det {Gah + Fab), (2.7) 
where Fab is the Abelian field strength tensoi[^ and Gab is the induced metric defined 

The Wess-Zumino part of the D7 brane action, Swz^ is given by 

Swz = Nffij I V[Y, A e^, (2.9) 

where 'P[C(p)] denotes the pull back of the background p-form field C(p) and = 1 + F -|- 
1/2F A F + . . . . The integral in (2.9) singles out the correct p-form, i.e. in our case a 8-form 
with legs along the world volume coordinates. 

For the embedding considered in (2.4) with all field strength tensors Fab = 0, the WZ 
action does not contribute and the DBI action Sdbi = —NffiTYo\{S^) f dt^xdr C reads 



~ a^a a^b (2- 



C = ^-det{Gab) = 1 + {drXsy + ^^^^('^-a^s)^. (2.10) 

A holographic topological insulator is given by a solution of the equations of motion fol- 
lowing from the Lagrangian (2.10) with the following boundary condition in the asymptotic 
AdS region 

lim Xs{r, xs) = Mq^ + O ( \] . (2.11) 



^In order to restore the correct 27ra' factors in the DBI and WZ action we have to rescale F — > 2Tra'F 
We do not distinguish between Xa and 
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Since Xs is interpreted as the mass of the flavor degrees of freedom is can be identifled 
with the position dependent mass m{x) for the fermions (and their superpartners) in (1.1). 
The equations of motion are quite comphcated and an exact solution with these boundary 
conditions is not known at present since. In [1] numerical solutions of the equations of motion 
were found using the heat method. 



3 Ansatz for a supersymmetric topological insulator 

In this section we show that the position dependent mass of the fleld theoretic topological 
insulator model presented in the previous section breaks all supersymmetries. However adding 
a localized counterterm can restore some of the supersymmetries. This argument motivates 
the modification of the holographic model by introducing a nontrivial gauge field on the probe 
D7 brane. 



3.1 Field theory analysis 

The low energy theory of the D3/D7 brane intersection considered in the previous section 
is given by = 4 supersymmetric Yang- Mills theory with gauge group SU{Nc) coupled to 
Nf M = 2 supersymmetric hypermultiplets. The A/" = 4 supersymmetric vector multiplet - 
describing the low energy theory of strings ending on the Nc D3-branes - can be decomposed 
into one M = 1 vector multiplet and three M = 1 chiral multiplets called $i, $2 and $3. The 
lowest components of the chiral multiplets, i.e. the scalar fields 0i, 02 and 03 can be identified 
with the six coordinates X4,X5, ... 5X9 transverse to the D3-branes: 0i ~ X4 + iX^,(f)2 ~ 
Xq + iXj and 03 ~ Xg + iXg. 

The low energy theory of the strings connecting the D3-branes and the Nf D7 branes 
are given in terms of Nf M = 2 supersymmetric hypermultiplets. For simplicity, we restrict 
ourselves to the case Nf = 1 in the following discussion. An M = 2 supersymmetric hyper- 
multiplet contains four fermionic and four bosonic degrees of freedom. It can be decomposed 
into an A/" = 1 chiral multiplet Q and an A/" = 1 anti-chiral multiplet Q\ or equivalently into 
two A^ = 1 chiral multiplets Q and Q. The scalar fields of these chiral multiplets are denoted 
by q and g, respectively, while the fermions are denoted by ■0 and ip. 

In the A^ = 1 language, the interactions of the vector and hypermultiplets are specified by 
the superpotential W 

W = e'^hi ($i$,$fc) + Q(m + <I>3)Q. (3.1) 

In the following only the superpotential term involving the mass of the quarks, m, will be 
important. Although not obvious in the A^ = 1 superspace language, the field theory is 
M = 2 supersymmetric and has a SU{2)-ji R symmetry as well as a global SU {2)l symmetry. 
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Note that both symmetries are also present in the brane intersection giving rise to the 5*0(4) = 
SU{2) L X SU{2)ti rotational invariance in the (X4, X5, Xg, X7) space. Note that the pair (g, g*) 
transforms as (0,1/2) under SU{2)l x SU{2)ti while the fermions [ipjip"^) are singlets under 
SU{2)l X SU{2)n. 

The realization of the topological insulator makes the mass of the fermions dependent 
on one of the spatial coordinates, where we choose 0:3 without loss of generality. This can 



be achieved by replacing m by m(x3) in the superpotential (3.1), introducing a position 



dependent masses for both the fermions and the scalar superpartners. It was shown [30] that 
position dependent couplings in the superpotential break all of the supersymmetries since the 
supersymmetric variation of the Lagrangian is no longer a total derivative. 

In [30] is was also shown that half of the original supersymmmetries can be restored by 
adding a counter-term to the action of the form (we give a review of the argument in appendix 
Bl) 

C^C + AC, (3.2) 



where AC is given by 



AC 



-2 Im 
. dm 



dm 6W \ 
dxs 5m{x-i) J 

^ 6W \ _ f m Y 

\bm[x-i)) ybmix'i)) 



(3.3) 



In the last step of equation (3.3) we assumed that m(x3) is real. Note that if the position 
dependence is of the form (2.11), i.e. m(x3) = M0X3/IX3I, the counterterm ( |3.3 ) will be 
delta-function localized at the interface location X3 = 0. For a real mass ^(0:3) and for the 



superpotential (3.1) the counterterm ( |3.3[ ) becomes 

AC 



(9X3 



(3.4) 



Let us compose the scalar fields q and q into a vector q (whose components are denoted by 
with m = 1, 2) 

*=(|.). (3.5) 
We can write the part of the counterterm quadratic in the scalar fields q and q as 



AC 



.dm ^ 

m-q q) 

0x3 

■ dm ^1 2 

g^a g + . . . , 

0x3 



(3.6) 
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where = q*'^ = {q*, q) and are the Pauh matrices. The terms represented by . . . in 
equation (3.6) are quadratic in the fermions if) and ip. 

Recall that field theory has SU{2) ^ x SU{2)n symmetry for a massive hypermultiplet and q 
transforms under (0, 1/2). Hence the counterterm ( 3.6[ ) has the following important properties: 
It has scaling dimension A = 2 and transforms as (0, 1) under the SU{2)l x SU{2)-ji symmetry. 

The two components of the four JV = 1 supersymmetry which is preserved by the ad- 
dition of the counterterm is given by (B.6). Since our field theory model is really J\f = 2 
supersymmetric the counterterm preserves four supersymmetries. 



3.2 Modified holographic ansatz 

The field theory analysis showed that in order to restore some supersymmetry a counterterm 
has to be introduced. In this section we present an ansatz for the holographic realization of 



such a counterterm. We need to identify a world volume field which can reproduce (3.3) at 
the AdS boundary. The counterterm corresponds to an operator of dimension A = 2 and it 
transforms under the SU{2)l x SU (2)7^ symmetry in the spin (0, 1) representation. In [31] the 
excitation spectrum on the probe D7 brane was analyzed in detail. It was shown an operator 
with these properties is dual to a KK excitation of the gauge field called (pj . In the following 
we give an explicit realization of that particular gauge field excitation. 

As we shall show the counterterm can be realized by turning on a world volume gauge field 
A which gives rise to a non- vanishing field strength tensor with components Fab = daAh — dbAa- 
The ansatz for the gauge field we choose is given by 



A 



\Ar J 



( X, \ 

— X4 
— X7 

V ^6 / 



2/5 

2/6 

\y7 J 



(3.7) 



In the second equality we have defined a second set of coordinates yi will be useful in the 
calculation of the BPS projectors later on. Note that the dependence of A on the coordinates 
Xa, a = 4,5,6,7 breaks the 50(4) = SU{2)l x SU{2)r down to SU{2)l. 



The relevant parts of the D7 brane action are given by 



S, 



D7 



-Nffxj j (fx ^-Aei{gab + F,,) + ]^Nf^y j P[C(4)] ^F ^F, 



where C(4) is given by 



+ X| + XlYdt A dxi A dx2 A dx^i + 



(3.9) 
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The . . . represent terms in C(4) which do not have legs along t, xi,X2, but have to be present 
in order to guarantee the self-duality condition F(5) = *-F(5) for the field strength F/^\ = dCu). 



For the gauge field ansatz (3.7) the Wess Zumino term of the D7 brane action reads 



Swz = -2^nYo\{S^) j d^xdr (r^ + Xlffdrf. 



(3.10) 



The DBl term of the D7 brane is not illuminating and will be evaluated later in (4.9) and in 



(4.10). 



Let us first study study the effects of the gauge field ansatz (3.7) and in particular determine 



the conformal dimension of the dual operator. Therefore we study the linearized equation of 



the gauge field for a D7 brane located at Xg = L,Xq = 0. With the ansatz (3.7) the term 



quadratic in / of the action for the Born-lnfeld part is given by (converting to polar coordinates 
which introduces a term in the integration measure) 



SDBi,,uad = -Nfii,^fo\{S') Jd'xdr['^{L' + r')'f + K{d,f)' + K{L' + ry{drfY], (3.11) 



while the WZ term (3.10) gives 



5, 



WZ 



-Nffij Yo\{S^) J d^xdr 2{L^ + r^ffdrf 
-Nfi^j Yo\{S^) J d^xdr {L^ + r^yOrif) 
4:Nffxjvo\{S^) J d'^xdrr{L^ + r'^)f. 



The equation of motion of / for the action Sdbi, quad + Swz becomes 



-d,.{r{r^ + L^fdrf) + {d,ff-A 



2\2 



-/ + 8(r2 + L2)/ = 0. 



Near the AdS boundary the asymptotic behavior of f is given by 

hm /(r) = c+(x)4 + c4x)^^ + ((r'^) . 



(3.12) 



(3.13) 



(3.14) 



By the standard holographic dictionary, turning on c+(x) corresponds to turning on a source 
for a operator of dimension A = 2. 



We conclude that the gauge field ansatz (3.7) transforms as (0, 1) under SU{2)l x SU{2)ii 



and has conformal dimension A = 2. Since it is the only type of fiuctuation with these 



properties, c+(x) as defined by equation (3.14) sources the counterterm (3.6). 
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3.3 Equations of motion 



The gauge field (3.7) is clearly not the most general form and one has to check that our 
ansatz is consistent. This means that the equations of motion for the gauge fields allow for 
all other components to be set to zero consistently. This is equivalent to the statement that 
all components of 

J d'x[d,(^^/-det{G + F){G + + (^d,{r^ + XD^y'^'^^A.^A" = (3.15) 

can be reduced to a single partial differential equation for /(r, x^). One can indeed show that 



the 6 = 3 component of the bracket in (3.15) is identically satisfied, whereas the a = 4, 5, 6, 7 
components are given by 

J^(^f,Xs,df,dXs,d^f,d^Xs)y'' = 0, (3.16) 



where y"' is defined in (3.7). The function J-" is a scalar function which is invariant under 
5*0(4) rotations in the 0:4, X5, Xg, space. Furthermore it depends only on /, Xg and their 
derivatives with respect to r and X3. If one expands the variation 6 A"" in terms of spherical 
harmonics, all of them vanish but the one proportional to ( |3.7[ ). This implies that there is no 
linear coupling of (3.7) to other spherical harmonics in the actioij^and hence the truncation 



is consistent. 

Both the equation of motion for Xg and / are highly nonlinear coupled second order partial 
differential equations whose form is very complicated and not very illuminating. We will not 



display their explicit form in this paper but they can be easily derived from the action (3.8). 



4 /^-symmetry for the D7 brane probe 

The K-symmetry transformation on the world volume of the D7 brane takes the form 

6e = {l + T)K (4.1) 

and is responsible for gauging away half the spinor degrees of freedom ensuring space time 
supersymmetry. The condition for an unbroken supersymmetry for the D7 brane probe is 
given by 

(l-r)e = 0, (4.2) 
where e are Killing spinors generating the 32 supersymmetries in the AdS^ x background 



[32] . It was shown in [33] (see also [31]) that for the AdS^ x space parameterized by (2.2) 



16 of the supersymmetries the AdS Killing spinors can be written in terms of constant spinor 



^Such a statement can also be derived from purely group theoretical arguments. 
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e = (r^ + p')'^6o, (4.3) 
where e is a doublet of ten dimensional Majorana Weyl spinors which satisfy 

ia2 ® roi23 Co = eo- (4.4) 



The spinor (4.3) generates the sixteen Poincare supersymmetries. In addition there are sixteen 
superconformal supersymmetries. For an embedding with nonzero Xg conformal invariance is 
broken and we only need to consider the Poincare supersymmetries. Since the r dependence 
of (4.3) is an overall factor, we can drop the overall factor in (4.3) and work with a constant 
spinor instead. 

The form of the /t-symmetry projector for the D-brane probe is given by [32]. 

v/-det(G + F)^2"n!^ ''^^^''^^^ ' ^^"'^ 

where 

jin) ^ (_l)n(^3)n ^ 701234567- (4.6) 

The gamma matrices 7^ are the pull-backs of the tangent space gamma matrices on the world 
volume and given by 

7. = E^d.X'^TA. (4.7) 

The explicit expressions for 7^ can be found in appendix [C] For the case of the D7 brane and 
a gauge field ansatz (3.7) the expansion of the projector in the field strength terminates at 
quadratic order. The /t-symmetry projector (4.5) can be expressed in terms of three prices 
Vn=k, /c = 0, 1, 2 with zero, one and two powers of the gauge fields respectively. 



r — r„=o + Tn=i + r 
1 



^-det(G + F) 



n=2 

icr2 ® 701234567 " (7^"^3a + ^l"'' Fab) CT i ® 701234567 
+ (^J^^'^FsaF,, + ^r'^'FabF,,] la^ ® 701234567 



(4.^ 



The determinant — det(G' + F) contains the metric determinant of the three sphere which can 
be extracted, where we define 

-det(G + F) = D (4.9) 

and D is given by 



D 



-2dsfdrfdrXsdsXs + (t' + {d,ff + (r^ + Xl)\dr.ff)Y (4.10) 
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4.1 Projection conditions on the spinors 



In order to simplify the form of the K-symmetry projector and solve the conditions for the 
existence of unbroken super symmetries we will impose three projection conditions on the 
infinitesimal supersymmetry transformation parameter e. 



i(^2 ® roi23e 
h ® r4567e 

cTs (g) r38r45e 



(4.11) 
(4.12) 
(4.13) 



Note that the first projection (4.11) was already imposed by singling out the Poincare super- 
symmetries (4.4). 

Since the three matrices on the left hand side all commute, the three projectors are com- 
patible and reduce the 32 supersymmetries of the AdS^ x background to 4. The expres- 
sions for r„=o,i,2 are evaluated in detail in appendix D, where the final form of the projector 
r = r„=o+rn=i+rn=2 can be found in (D.18). We display this result again for the convenience 
of the reader: 



Te 



1 I ^r^S r 



xlffd,f-d,x,[ 



2/ , drf 



le 



+ 



+ 



ir' + XI) 



-d.X. 



+ X|)2 r 
-j= i OrXs >l2®x Tdse 



^ + — ) + ^ {drfdsXs - dsfdrXs) } h ® T,s e 



1 Udsf 



Dr 



-drXsh2(^x''Td8e. 



(4.14) 



It is remarkable after employing the conditions (4.11 )-(4.13 ) the complete projector can be 
written as the sum of four terms involving only four linearly independent combinations of 
gamma matrices. The conditions (4.11)-(4.13) reduce the Poincare supersymmetry from 32 to 
4. Consequently the ^-symmetry projector (4.2) should not reduce the number of unbroken 
supersymmetries further. 

It follows that r has to be equal to the unit operator. This implies that the term in (4.14) 
proportional to I2 ® 1 has to be equal to one, or equivalently 



1 + ^3/ 



+ XirfdJ-dsXs{% + ^) = ^^. 



(4.15) 



where D is given by (4.10). Furthermore, the linearly independent terms proportional to 
® x'^Tds and I2 ^ x'^T^s have to vanish separately. Note that the vanishing of 



1 



r 



38) 



terms proportional to I2 ® x'^Tds and I2 ® x'^Tds give the same condition, hence we have two 
additional BPS equations 



12 



^3/ 



drXs 



'2/ , drf, , 2/ 



- + — ) + ^ {drfdsXs - dsfdrXs) 



(4.16) 
(4.17) 



A lengthy calculation shows that if (4.16) and (4.17) are satisfied then both (4.15), as well 



as the equations of motion, are automatically satisfied. Hence finding a D7 brane embedding 
which preserves four supersymmetries boils down to finding solutions two BPS equations (4.16) 
and (|4l7|. 



5 Solutions of the BPS equations 



The main result from the last section was the the there are four unbroken supersymmetries 



if the BPS equations (4.16) and (4.17) are satisfied. The two equations are considerably 

They are 
In the 



simpler than the equations of motion, since they involve only first derivatives 
however still nonlinear partial differential equations and hence not easy to solve 
following we will present some solutions to the BPS equations using a variety of methods. 
The solution corresponding to a supersymmetric holographic topological insulator will be 
obtained numerically. 



5.1 A solution using scaling symmetry 



The two BPS equations (4.16) and (4.17) are invariant under the following scaling transfor- 



mation (which generalizes the scaling transformation of [T]). 



(5.1) 



One can use the scaling symmetry to obtain an ansatz for the fields which only depend on the 
scaling invariant combination y = x^r 



Xs = rh{y), f = f{y). 



(5.2) 



The BPS equations become 



h' 



yh' + h-f = 0, 
+ 2f + {y + 2fh)f = 0. 



(5.3) 
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Figure 1: (a) numerical solution for f{y), (b) Plot of f{r,x^), (c) Plot of X^{r,x^) 



The first equation can be integrated and determines h 

y y 



(5.4) 



The second equation of (5.3) then takes the following form 



(y2 + 2/(ci + /))(l+(^ 



2 • 



(5.5) 



While this ordinary differential equation does not seem to have an explicit solution in terms 
of known functions, it can be easily integrated numerically. Note that if one imposes the initial 
condition /(j/)|j/=o 7^ the solution indeed satisfies 9^X81^=0 = which is equivalent to the 
condition that the 3-sphere closes off smoothly at r = 0. In figure [T] we display an example 
for ci = and f{y)\y=Q = 2. 

Note that this solution does not correspond to a topological insulator. This follows from 
the fact that lima;3_!.-|-oo = which means that the fermion mass cannot behave as in (2.11 ). 



5.2 Perturbation theory 



In this section we will give a perturbative solutions of the BPS equations (4.16) and (4.17). 



First, we redefine f = rf such that the BPS equation (4.16) becomes 

d^f = drXs 



while the second BPS equation (4.17) reads 

dsXs 



drf+-f + 



+ % ( rd^fdrXs - d,Xs{f + rdj) ) = 0. 



(5.6) 



(5.7) 
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Equation (5.6) can be solved by 

/ = drh + d^h, Xg = dsh — drh, 



(5.8) 



where h can is a harmonic function satisfying 

{dl + dl) h = 0. (5.9) 

We should consider h as a fluctuation on top of the embedding. Since we are only interested 
in the embedding, we set h = from now on and plug equation (5.9) into the second BPS 
equation (5.7). This leads to a complicated partial differential equation for the function h 



dlh 



+ 



3(9,/i 2dlh {drhf 2drh{drd^hf 2dlh d^h drh 



(r2 + (93/i)2)2 r 



+ 



+ 



+ (9,2/1 = 0. (5.10) 



Unfortunately, we were not able to solve equation (5.10) analytically. Following fTj we take 
the ansatz for the perturbation expansion for h as follows 



h = e-hi{rx^) + e'—h2,{rx^) + 



n=l 



2n-l 



The order e contribution to (5.10) is given by, renaming rxs = y 

h{y) - vKiy) - (i + y')K{y) = 0, 

which has the most general solution 

/ii(y) = Mo(l + 1/2)1/2 + cii/. 



(5.11) 



(5.12) 



(5.13) 



In the following we set ci = since it corresponds to a constant mass term. In this case the 
order contribution to (5.10) has the form 



31/ 



1 + 1/2 2 ' l^y2 

This equation can be integrated and one gets 



(5.14) 



l+6i/2+V)+C2,ii/+C2,2(3l/ arcsinh(y) + v/l + 1/2(1/2-2)), (5.15) 



where the second and third terms are solutions to the homogeneous part of the (5.14). We 
note that the term of order is dominant near the boundary, i.e. as r — )■ oo one obtains from 

+ 0{-^ \ =MoP^ + 0{ ^ 1 . (5.16) 



(|5l3| and 



lim Xg = lim Mq 



^/lT~(x^ 



F3 
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Figure 2: (a) numerical solution for X8(r, X3), (b) Plot of f{r,x^) 

for large r. 



It is in principle possible to find the higher orders of the perturbation series (5.11), which 



are all sub-leading in powers of r. However the perturbative solution shares the same problem 
as [1] . For the leading term in e one has 



lim drX^ 



M0X3 + 



(5.17) 



One can easily convince oneself that the non vanishing of drX^ at r = is a general feature 



of the perturbation series (stemming from the scaling symmetry of section 5.1). However 
drX^ has to vanish at r = for the sphere to close off smoothly. Consequently, while the 
perturbation expansion has a reasonable behavior at r — )■ 00, the behavior as r — )■ does not 
correspond to a smooth solution. 



5.3 Numerical solution of the BPS equations 



In this section we solve the BPS equation (5.10) numerically by using a heat method. This 



method was also used to obtain solutions for the non-supersymmetric setup of [T]. 

In particular we make /i(r, X3) dependent on an imaginary time r. The derivative with 



respect to r is given by the left-hand side of equation (5.10) multiplied by r 



dh{r,X3,T) 
Or 



(r2 + {d^hyy 

+2r drKdrd^hf 



+ r'dih + 3r drh + 2dih {drhy 



2r dlh dlh drh. 



(5.1^ 
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After a sufficiently long time r the field configuration approaches a solution of the right 
hand side equation. We solved this equation for a spatial region X3 = [—Xmax,Xmax] and 
r e [0, Tmaa;]- As a good initial function h at time r = we use the first order solution (5.13) 
(with Ci = 0) which we can trust for large r. At the boundaries of the spatial regions, we impose 
Neumann boundary conditions for h at r = and Dirichlet boundary conditions (given by the 
first order solution). The numerical solution of the equations was obtained using Mathematica 
and unlike the perturbative solution the numerical solution has the correct behavior both in 
the r = region and in the large r region. 



6 Discussion 

Making coupling constants dependent on space time coordinates generically breaks all the 
supersymmetries of a supersymmetric field theory. Examples are Janus solutions [32] and 
the holographic topological insulator solutions constructed in [1] . Turning on a world volume 
gauge field allowed us to preserve four of the sixteen supersymmetries of the probe D7 brane. 
On the field theory side this corresponds to turning on localized counterterms as expected 
from the weak coupling analysis. A supersymmetric generalization of the topological insulator 
might be interesting in its own right from a theoretical perspective. 

The BPS equations we have found are considerably simpler than the equations of motion, 
since they are first order instead of second order partial differential equations. Utilizing the 
scaling symmetry we have been able to generate exact smooth solutions. These solutions do 
however not have the boundary behavior which describes topological insulators. We have been 
unable so far to obtain closed form solutions corresponding to holographic topological insula- 
tors. Instead we had to resort to other methods such as perturbation theory and numerical 
methods. A closed form solution would be very useful as it is possible to study fluctuations 
about the solution and calculate correlation functions and transport quantities. 

It is an interesting open problem to determine whether the difficulty of obtaining closed 
form solutions is a fundamental problem due to the small amount of preserved supersymmetry 
or simply due to the lack of ingenuity on the authors' sid^ 

Even with the numerical solutions it might be feasible to extract observable quantities, like 
transport quantities from the numerical solutions. We leave such questions for future work. 



'^Note that in the case of supersymmetric Janus solutions [M] the original form of the BPS equations looked 
too complicated to be solved. Only a sequence of clever variable changes made a solution possible. 
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A Conventions 

We will briefly review our conventions for the IIB supergravity and ten dimensional gamma 
matrices which follow |32]- The supersymmetry transformation parameters e of type IIB 
supergravity can be written as a doublet 

-(::)■ (-^) 

The flat space gamma matrices Fj are 32 x 32 matrices and the spinors e^, i = 1, 2 are Majorana- 
Weyl spinors 

rnei = ei, z = l,2. (A.2) 
The matrices in the K-symmetry projectors are of the form 

CTj ® ri^j2...i„, (A. 3) 



where Ti^^i^.i^ = Ti-^Ti^ . . . Fj^ and the Pauli matrices act on the doublet (A.l) 



B Supersymmetry and counter-terms 

In this appendix we review and adapt the results given in for the field theory description 
of the topological insulator discussed in the body of the paper. 

The theory is most conveniently formulated in terms in = 1 chiral multiplets, where 
only the A^ = 1 supersymmetry is manifest. We consider n chiral multiplets denotes with 
the field content. 

$V (f)\ F\ t = l,2,--- ,n, (B.l) 

where 0* is a complex scalar, ip^ is a Majorana spinor and is an auxiliary field. In compo- 
nents the Lagrangian density is given by 

i -■ ■ ■ i dW -■ dW ■ 
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where P+ = |(1 + 7^). The supersymmetry transformations are given by 

6F' = -iV2e^^'d^P+tl)\ (B.3) 

Here e is a Majorana spinor parameter. W is the superpotential which parameterizes the 
self interactions of the chiral superfields. It was shown [30] that couphngs in the superpo- 



tential which position dependent break all of the supersymmetries (B.3). Assuming that the 



superpotential depends on one coupling constant g{x^)^ the supersymmetry variation of the 



Lagrangian (B.2) is given by 

is not a total derivative. Consequently, all supersymmetries are broken by the position de- 
pendent mass. Furthermore it was shown in [SU] that some - but not all - supersymmmetries 
can be restored by adding a counter-term to the action of the form 

dg \fdW\ fdWy] , 

The supersymmetry which is preserved is given by 

ne = e, n=^(l + z7V)- (B.6) 



For the theory discussed in section 331 the coupling g is identified with the mass m, the chiral 



superfields $1, $2 are identified with Q,Q the superpotantial is given by 

W = m{x3)QQ (B.7) 

C D7 brane embedding 

In this appendix we present the details of the D7 brane embedding for completeness. The 
AdS^ X background metric is given by 

ds^ = gf^udx^dx" = (r^ + p'^)r]ijdxidxj H — ^ ^ j dxl + dX^ + dXg j , (C.l) 

^ \a=4 J 

where 

= xl + xl + xl + xl p2 _ xl + Xl (C.2) 
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Note that we have introduced two different set of indices: i,j G {0, 1, 2, 3} and a G {4, 5, 6, 7}. 
The corresponding coordinates Xi and Xa parametrize the field theory coordinates and the 
internal coordinates, respectively. The D7 brane world volume coordinates ( are identified 
with Xi{i = 0, 1,2,3) and Xa{a = 4,5,6,7). In order to embed the D7 brane into AdS^ x 
we have to specify Xg and Xq as a function of the world volume coordinates. In this paper 
we restrict ourselves to 

Xs = X8{x3,r), 

Xg = 0. (C.3) 

The induced metric on the world volume of the brane, with components Gab, is defined by 

Gab = daX^dbX'g,, (C.4) 

and takes the following form 

Gij = {r^ + Xi)rjij, 

G33 = ir' + Xi) + -^^^d,Xsd,Xs, 



G. 



}_ Xa 

+ Xl r 
1 



^3a — ..9 , v2 ~7^2,X^drX^, 



ab 



r2 + X| 



[Sab+idrXsr"^), (C.5) 



where i, j G {0, 1, 2} and a,b E {4, 5, 6, 7}. The non- vanishing components of the field strength 
are, for a = 4, 5, 6, 7 

F, = dsf{r,xs) 
fir, X3) 

= -'^-h[l^)ixl^xl), 

1 f fir,xs) \ , , 

F47 = -Or I ^ I {X4X6 - 2:5X7), 

F56 = -Or I ^ 1 {X4Xe - X5X7), 

F57 = ( ^^^^2^"^ ^ {X4X7 + X5X6), 

2f{r,X3) I f f{r,X3)\ 2 . 2^ 

= —2 + -dr\^—^j{Xe + X,). (C.6) 
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Note that F^q = — -F57 and F47 = F^q. The pull back of the gamma matrices on the world 



volume (4.7) are given by 



73 
7a 



r2 + x|r„ 



0,1,2, 



H + x«2r3(i + -^r 



1 



r2 + X, 

H XaTs) 



2^ 38 



4,5,6,7. 



(C.7) 



D Calculation of the At- symmetry projector 

In this appendix we will present the detailed calculation of r„=0) ^n=i and r„=2- We use the 



three projection conditions (4.11 )-(4.13) as well as the properties of the world volume gauge 



field to bring the projector in a minimal form. 



D.l Calculation of r„=o 

Using the expressions for the pulled back gamma matrices one get 

1 



n=0 



V-det(G + F) 
1 



(X> 701234567 



ia2® roi234567 



1 + / 2 1"^ v2\ ^38^3-^8 + ^^(a;ar„)r8 



Employing the projectors ( |4.11[ ) and ( |4.12[ ) r„=o can be expressed as follows 



D 



r2 + Xi) 



dS 



e. 



D.2 Calculation of r„=i 

The contribution to F which is linear in the field strength is given by 



7'"i^3a + ^Y'Fab ] (Tl ® 701234567 



(D.i: 



(D.2) 



(D.3) 



where e"'"^'^ is the totally antisymmetric tensor in (4567) space normalized such that e 



4567 



1. 
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We can use the formulae for the pull back of the gamma matrices (C.7) to determine 



70123c(i 



(r^ + Xg)roi23cd + d3XsToi28cd + {r"^ + Xg)-^—^roi23{Xdrc8 - XcTdS), 



d Xg 

10l2bcd = ^012bcd H {XbTQi2cd^8 — 3:^012^1^8 + XdToi2bc^8) 

r 

In addition we need the following relation involving the field strength 

4/ 



-Vd^d e, 



(D.4) 



(D.5) 



where the coordinates yi were defined in (3.7). There are three additional identities involving 
the field strength, which we need in the following 

■1 



/ , 2a,/ 



— Fg? )e 



as well as 



and 



e'^'^^F^aXi^T^d^e = -53/(145 - T^^y^ 



e -Tsa-L bed e 



(D.6) 



(D.7) 



(D.8) 



where we used the projection condition (4.12) to simplify the expressions above. Employing 
these identities together (D.3) can be expressed as follows 

2 

= (r" + X^)- — -^Vria^ (X) rm9sr^s + - ^ 



rn=i e 



D 



+ Xl)^^y,ai ® Toiasr^s + ^(^' + ^s) (4 + ® roi23(r45 - Tst) 



=53X3(4 + ® roi28(r45 - Tst) + ^^as/Ai ® roi2rf 
53/0-1 (g) roi28(r45 - Te?) 



D 

1 la^Xg 



v^2 r 



(D.9) 



This can be simplified by using (|4.11|), 
2 



8M^2 ' 2r 



+ X^ 



^3^ydrd8 + ^ir' + x?^^^ ■ 



D 



0-3 ® (r45 - Ter) 



1 1 



^53X3(4 + 1^)^3 ® r38(r45 - r67) - ^^53/1/^3 ® r3, 

'D V^" 2r / V-D ^ 



•93/^3 O r38(r45 - Tst) 



1 1(9,X8 
D2 r 



(D.IO) 
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The dependence on the can be removed by employing (4.13) from which the following 
identities can be derived 



CTg (g) (145 - r67)e 

0-3 (g)r38(r45 - r67)e 



Then the r„=i projector becomes 



D 



8'' ^ ^2 



+x'^ I2 ® Tdse, 
-2 la^Tsse, 
2 I2 ® le. 



^,3-— (r' + X|^ 



r / \/Df 



+— = dsf I2 ® 1 



e. 



(D.ii: 



■' +^)l2®r38 



(D.12) 



D.3 Calculation of r„=2 

The part of the projector which is quadratic in the field strength is given by 



1 /I 



( ^7'"''i^3aFfec + T:r'"'''FabF,dUa2 ® 701234567 



D\2 

^ : (le^''"'F3aFi„ia2 ® 7oi2d + le''''"'FabF,d^(T2 ® 70123 



Expressing the world volume gamma matrices in terms of fiat space ones gives 



(D.13) 



70123 = (r^ + X8^)%i23 + (r' + ^8)^3^8roi28, 
7012d = {r^ + -^8 ) (roi2d + -^XdToi28 



(D.14) 



and using the following relations which follow from the explicit ansatz for the gauge field (3.7) 

2 



8 



abed 77' 77' 
e i'ab-fcd 



1 



abed p p.V, 
^ 3a-' be'- a 



+ 4/^3/ 



le^'^'FsaFbcXd = +^fdsf, 



(D.15) 
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we obtain 
rn=2e 



= + XlYfdrf la^ ® roi23 - 4(^' + XDfdrfd^X, lO^ ® TqisS 

(D.16) 



Using the projector (4.11) the expression can be simphfied further to give 
1 f 2 



+ X|)2/9./ I2 ® 1 + ^(r^ + XDfdrfd^X^ I2 ® r 
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D y 

+ -(r^ + XD/^s/ x'^ I2 ® Trfs - ^(r'^ + XD/^s/ a.Xg I2 ® Tag ^ e. 



(D.17) 



D.4 Final form of the F projector 



The expressions for r„=o,i,2 in are evaluated in append D The final result for the projector 



r = r„=o + ^n=i + ^n=2 from (D.2), (D.12) and (D.17) one gets the final form of the V 
projector. 

re = ill + ^a3/-4(r2 + X|m./-93X8(^ + ^)|l2®le 



r 



+ 
+ 



(r^ + XI) 


r -^3X8 




I (r2 + X|)2 



-{% + —) + % {drfd3Xs - dsfdrXs) } h ® T^S £ 



2{r^ + Xl)f ^d,f 



drXAU®x''Vd.e + 



1 ir^s/ 



(D.18) 
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